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ABSTRACT 

We present optical spectroscopic data of the inner circumstellar ring around SN 1987A from the 
Anglo-Australian Telescope (AAT) and the Very Large Telescope (VLT) between ~1400 and ~5000 
days post-explosion. We also assembled the available optical and near-infrared line fluxes from the 
literature between ~300 and ^2000 days. These line light curves were fitted with a photoionization 
model to determine the density structure and the elemental abundances for the inner ring. We found 
densities ranging from lxlO 3 to 3xl0 4 atoms cm -3 and a total mass of the ionized gas of ~5.8 x 
10~ 2 M© within the inner ring. Abundances inferred from the optical and near-infrared data were also 
complemented with estimates of Lundqvist & Fransson (1996) based on ultraviolet lines. This way we 
found an He/H-ratio (by number of atoms) of 0.17 ± 0.06 which is roughly 30% lower than previously 
estimated and twice the solar and the Large Magellanic Cloud (LMC) value. We found an N/O-ratio 
of 1.5 ± 0.7, and the total (C+N+0)/(H+He) abundance about 1.6 times its LMC value or roughly 
0.6 times the most recent solar value. An iron abundance of 0.20 ±0.11 times solar was found which 
is within the range of the estimates for the LMC. We also present late time (~5000 - 7500 days) line 
light curves of [O III], [Ne III], [Ne IV], [Ar III], [Ar IV], and [Fe VII] from observations with the VLT. 
We compared these with model fluxes and found that an additional 10 2 atoms cm~ 3 component was 
required to explain the data of the highest ionization lines. Such low density gas is expected in the H 
Il-region interior to the inner ring which likely extends also to larger radii at higher latitudes (out of 
the ring plane). At epochs later than ^5000 days our models underproduce the emission of most of 
these lines as expected due to the contribution from the interaction of the supernova ejecta with the 
ring. 

Subject headings: circumstellar matter — supernovae: individual (SN 1987A) 



1. INTRODUCTION 

Supernova (SN) 1987A, with its surrounding triple ring 
structure is one of the best known astronomical images 
from the 1990's. However, the origin of these spectacular 
rings is still somewhat uncertain. It is believed that the 
material in these rings was ejected by the SN progenitor 
star ^20,000 years ago. Also, the elemental abundances 
derived for the inner ring {e.g., Lundqvist & Fransson 
1996 hereafter LF96) strongly support a stellar origin for 
the ring material. Therefore, studies of the ring structure 
can also provide additional information on the progeni- 
tor star of SN 1987A which has been shown by detailed 
pre-explosion observations to have been initially a 14- 
20 Mq star and a blue supergiant of B3I spectral type 
at the time of the explosion (for review and references 
see Smartt et al. 2009). Recently, binary merger models 
(e.g., Morris & Podsiadlowski 2007) have proposed to ex- 
plain the blue supergiant nature of the SN progenitor as 
well as the formation of the triple ring system. A rapidly 
rotating single star progenitor has been suggested as an 
alternative origin for the ring structure (e.g., Chita et al. 
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2008). 

The gas in the ring structure was photoionized by the 
extreme UV and soft X-ray SN flash with the first nar- 
row optical and UV emission lines detected a few months 
after the explosion (Wampler et al. 1988; Fransson et al. 
1989). The gas then cooled and recombined until the 
start of the collision between the expanding SN ejecta 
and the material of the inner ring. The first evidence 
that the ejecta/ring impact had started came from ra- 
dio (Staveley-Smith et al. 1992) and X-ray emission 
(Gorenstein, Hughes & Tucker 1994) in mid-1990. The 
first signs at optical wavelengths were detected in Hub- 
ble Space Telescope (HST) images from 1995 (Lawrence 
et al. 2000). Recent high resolution X-ray images (Ng et 
al. 2009) from 2008 indicate that the SN blast wave has 
already overtaken the inner ring. The UV-optical emis- 
sion line fluxes from the inner ring have been modeled 
by LF96 up to ^2000 days from the explosion. They 
found gas densities ranging from 6 x 10 3 cm~ 3 to 3.3 
x 10 4 cm -3 and derived the following relative (ratio by 
number of atoms) abundances: He/H = 0.25 ± 0.05, 
N/C = 5.0 ± 2.0 and N/O = 1.1 ± 0.4 with an overall 
metal (i.e. C, N, and O) abundance of 0.30 ± 0.05 times 
solar. HST data up to ~ 3500 days were analyzed by 
Lundqvist & Sonneborn (1997) and showed evidence of 
densities in the ring down to ~2 x 10 3 cm~ 3 . Abundance 
estimates for metals in the SN 1987A CSM have recently 
also been made from X-ray observations (Zhekov et al. 
2006, Dewey et al. 2008, Heng et al. 2008, Zhekov et 
al. 2009, 2010) of emission originating from the interac- 
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tion between the expanding SN ejecta and the material 
in the inner ring. These observations yielded a similar 
N/O abundance as found by LF96, although the absolute 
C+N+O abundance is lower, and they also add abun- 
dance estimates of heavier metals. 

The main aim of the present study is to investigate 
the density structure and the elemental abundances of 
the inner ring using optical emission line fluxes from 
the entire era before the collision between the SN ejecta 
and the ring material started to dominate these fluxes. 
To do this, we obtained optical spectra of the SN 
1987A CSM with the Anglo- Australian Telescope (AAT) 
between — 1400 and —4300 days and with the Very 
Large Telescope (VLT) between -5000 and -7500 days 
from the explosion. We also collected all the available 
optical/near- IR absolute line fluxes of the circumstellar 
ring from the literature for epochs between —300 and 
—2000 days. This data-set covers epochs between 300 
days and —7500 days from the SN explosion i.e. the en- 
tire era from the time when the emission from the inner 
ring first became visible at optical and near-IR wave- 
lengths until and beyond the first signs of the collision 
between the SN ejecta and the inner ring were detected at 
these wavelengths. Preliminary studies of the line light 
curves between —300 and —4300 have already been pre- 
sented in Mattila (2002) and Mattila et al. (2003). The 
late time VLT observations have also been reported in 
Groningsson et al. (2008a, b). In Section 2 we describe 
the AAT and VLT observations, data reductions and the 
available line fluxes from the literature. In Section 3, 
the line light curves arc modeled with a photoionization 
model to estimate the density structure and elemental 
abundance of the inner ring. Finally, in Section 4 we 
discuss the results and in Section 5 present conclusions 
from this study. 

2. OBSERVATIONS 

2.1. AAT Spectroscopy 

Optical spectroscopic observations of SN 1987A were 
carried out using the Royal Greenwich Observatory 
(RGO) Spectrograph on the AAT between 1991 and 1998 
(for full details of the AAT observations see Table 1). 
Data were obtained at four different epochs: 1416, 1680, 
1991, and 4309 days post-explosion. At each epoch, the 
observations comprised a brief, wide slit integration, and 
a longer duration, narrow slit integration (see col. 6 in 
Table 1) with spectral resolutions of FWHM -300-600 
km s _1 and —150 km s _1 , respectively. The SN observa- 
tions spanned air masses —1.5 to —3.0, making accurate 
flux calibration quite challenging. We therefore carried 
out the observations with the slit position angle (PA, col 
5 in Table 1) set to be roughly the same as the line join- 
ing Stars 2 and 3 (see Fig. 1). This meant that Star 
2 lay well within the wide slit and so could be used to 
correct for the variable transmission of the atmosphere. 
However, an additional problem was that, since the slit 
PA was generally not at the parallactic angle, when the 
narrow slit was used, atmospheric refraction could intro- 
duce wavelength-dependent vignetting (e.g., Filippcnko 
1982) of the ring spectra. Moreover, the magnitude of 
this effect differed from that experienced by Star 2 since 
the latter lay nearer the edge of the slit. Consequently, 
the fluxing uncertainty introduced could be as large as 




Fig. 1.— HST WFPC2 F656N image of SN 1987A ring structure 
from 10th July 1997 (3790 days). The position angle (-60°) and 
width (5.3" corresponding to 1.3 pc at the distance of the LMC) 
of the wide slit used for the AAT observations on Dec. 1994 arc 
indicated. North is up and East is to the left. 



±40%. 

Data reduction was carried out using standard IRAF 
routines and for the 2D frames included the overscan 
correction, bias subtraction, flat fielding (to correct for 
fringing) and removal of cosmic rays. The background 
subtraction was performed using the IRAF BACKGROUND 
task in two steps, first to remove the emission of atmo- 
spheric lines and then the more complex Large Magel- 
lanic Cloud (LMC) background. The wavelength cali- 
bration was performed for the background subtracted 2D 
frames using CuAr arc lamp spectra observed together 
with the SN. An unweighted extraction was performed 
for the ring spectra because under the AAT seeing con- 
ditions the spectra of the two nearby stars were always 
blended with the spectrum of the ring making the use 
of any weighting algorithms unreliable. For the wide-slit 
observations, extraction apertures were large enough to 
include the flux from both the two nearby stars and the 
inner ring (see Fig. 1) i.e. all the ring flux was included. 
For the narrow slit observations smaller extraction aper- 
tures were used in order to maximise the S/N of the 
ring spectra and to avoid contamination by the strong 
stellar absorption lines in Stars 2 and 3 prominent at 
wavelengths shorter than —5000 A. The telluric correc- 
tion was performed using spectra of the low metallicity 
F5 star HD 26169 located close to SN 1987A on the sky. 
This star has only a very few stellar features apparent in 
its spectrum which were removed by interpolation. One 
of the spectrophotometric standards EG 21, LTT 1788, 
LTT 2415 or LTT 3218 (Hamuy et al. 1994) was observed 
with each of the SN observations and was used for the 
relative flux calibration of the SN spectra. Examples of 
the reduced narrow-slit spectra are shown in Fig. 2. 

The SN observations spanned airmass values between 
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~1.5 and ~3.0. At such high airmasses some of the ob- 
servations are likely to be affected by clouds if the night is 
not completely clear. Therefore, each wide slit spectrum 
was checked for the effects of attenuation by clouds, and 
a correction applied to the affected spectra. This was 
carried out by making use of the nearby Star 2 (for full 
details see Mattila 2002). The spectrum of Star 2 was 
extracted from each of the wide slit observations and 
its continuum was fitted with a low order polynomial. 
The continuum with the highest counts (which also had 
the expected shape for a stellar continuum) was selected 
to represent the unattenuatcd Star 2 continuum. To 
compensate also for possible uncertainties in the relative 
fluxing a correction curve was created for each wide slit 
spectrum by dividing the 'real' Star 2 continuum by the 
observed one. The extracted wide slit SN spectra were 
then multiplied by the corresponding correction curves. 
In the AAT spectra the narrow lines arising from the 
inner ring were interactively fitted with a Gaussian pro- 
file on top of the broad SN ejecta lines. As these lines 
have very different widths (e.g., see Fig. 2) the SN ejecta 
lines could be treated as a local background for the nar- 
row CSM lines and fitted with a low order polynomial 
on the two sides of each narrow line. We estimate that 
the uncertainties in the absolute fluxes measured from 
the wide slit spectra are ~±15%. For the narrow slit 
spectra there is also another source of uncertainty in the 
relative fluxes due to the differential refraction of light 
in the Earth's atmosphere. The flux calibration of the 
narrow slit spectra and estimation of the associated un- 
certainties are described in Sect. 2.2. 

2.2. Absolute Flux Calibration with Archival HST Data 

The AAT observations were obtained with the seeing 
ranging between ~1" and ^3" and so the inner ring spec- 
tra were always significantly contaminated by light from 
Stars 2 and 3 located 2.9" and 1.6", respectively, from 
the SN (see Fig. 1). In general, the CSM emission lines 
could be easily distinguished from the continuum orig- 
inating from the SN and the two stars. However, Star 
3 is a Be star showing strong Ha and H/3 emission to- 
gether with optical variability of ~0.5 magnitudes (Wang 
et al. 1992). Consequently, it could affect significantly 
the observed line fluxes from the inner ring, especially 
at the later epochs when the ring was fainter. To as- 
sess and correct for possible Star 3 contamination, we 
searched the HST archive for suitable optical spectro- 
scopic (STIS) and photometric (WFPC2) observations 
(see Table 2) covering the epochs of our AAT observa- 
tions. 

STIS spectra with a large enough slit aperture 
(52" x2") and suitable centering and orientation to in- 
clude Star 3 were selected from the HST archive at 
three epochs of observation. Pipeline reduced STIS data 
were acquired from the HST archive and tasks from the 
IRAF STSDAS package were used to combine the pipeline- 
reduced single exposures, rejecting for cosmic rays and 
hot pixels. The emission lines from the sky were sub- 
tracted using the IRAF BACKGROUND task with a 1st or 
2nd order polynomial fit to the sky lines. However, 
the subtraction of the background emission originating 
from spatially extended structures such as the outer rings 
was not possible and therefore, some faint nebular back- 
ground emission is visible in the extracted Star 3 spec- 



TABLE 1 
Log of the AAT spectroscopy 



Date 


Epoch Wavelength 


Exp. 




Slit width 


Seeing 




[days] 


[A] 


[sec] 


roi 


r//i 


r//i 
[ 1 


1991 


1416 


5030-6300 


2000 


-65 


6.7 


2.0 


Jan 9-10 




6175-7450 


1000 


-65 


6.7 


1.8 






7130-9900 


1000 


-65 


6.7 


2.0 






5030-6300 


5400 


-65 


2.0 


2.0 






6175-7450 


6400 


-65 


2.0 


2.0 






7130-9900 


6900 


-65 


2.0 


2.0 


1991 


1680 


6509-7790 


300 


-42 


6.7 


2.7 


Oct 1 




7600-8900 


300 


-42 


6.7 


2.7 






5400-6680 


300 


-42 


6.7 


2.6 






4314-5576 


300 


-42 


6.7 


2.6 






3220-4470 


300 


-42 


6.7 


2.4 






3220-4470 


3600 


-42 


2.0 


2.4 






4314-5576 


3600 


-42 


2.0 


2.6 






5400-6680 


3600 


-42 


2.0 


~2 






6509-7790 


3600 


-42 


2.0 


2.7 






7600-8900 


3600 


-42 


2.0 


^2 


1994 


2864 


3200-4780 


1000 


-60 


5.3 


1.1 


Dec 27-28 




5980-7600 


800, 100 


-60 


5.3 


1.1 






4600-6200 


1000 


-60 


5.3 


1.1 






7380-9000 


600 


-60 


5.3 


1.4 






3200-4780 


lUoUU 


fin 
-OU 


1.0 


l.l 






4600-6200 


10800 


-60 


1.5 


1.1 






5980-7600 


9000 


-60 


1.5 


1.1 






7380-9000 


7600 


-60 


1.5/1.7 


1.4 


1998 


4309 


3250-6370 


600 


-60 


10.0 


1.1 


Dec 10-11 




5180-8320 


600 


-60 


10.0 


1.1 






7850-11040 


4400 


-60 


10.0 


~3 






3250-6370 


7200 


-60 


2.0 


1.9 






5180-8320 


1200 


-60 


2.0 


1.1 






6000-7600 


8000 


-60 


2.0 


2.2 



TABLE 2 

HST OBSERVATIONS OF STAR 3 AND THE INNER RING 



Date 


Epoch 


Instr. 




Star 3 


Inner ring 




(days) 




Ha 


H/3 


m(B) 


F656N F502N 


1994 Feb 3 


2537 


WFPC2 


39 




15.97 


226 (244) 57 (70) 


1994 Sep 24 


2770 


WFPC2 






15.96 


48 (50) 


1996 Feb 6 


3270 


WFPC2 






15.81 


38 (41) 


1997 Apr 26 


3715 


STIS 


41 


12 






1997 Jul 10 


3790 


WFPC2 


45 




16.03 


209 (167) 33 (39) 


1997 Dec 9 


3942 


STIS 


46 


13 






1998 Feb 5 


4000 


WFPC2 


46 






175 (157) 


1999 Apr 21 


4431 


WFPC2 


52 






166 


2000 May 1 


4816 


STIS 


36 









The Ha and H/3 fluxes and £?-band magnitudes of Star 3 and 
F656N and F502N fluxes of the inner ring as obtained from the 
HST observations. The line fluxes for Ha (after correcting for 
Star 3) + [N II] A6548 and [O III] A5007 as obtained for the inner 
ring from the AAT observations and interpolated to the epochs of 
the HST observations are given for comparison in brackets after 
the F656N and F502N fluxes, respectively. All the fluxes are in 
the units of 10 -15 ergs -1 cm~ 2 . 

tra. The Star 3 spectra were extracted using the STSDAS 
X1D task with the recommended aperture size of 5 pix- 
els (0.25") selected to match the size of the stellar PSF. 
This task traces the center of the spectrum by cross cor- 
relation and automatically performs a 2-D wavelength 
calibration and flux calibration for the data. The only 
strong lines visible in the reduced Star 3 spectra are Ha 
and H/3 confirming the findings of Wang et al. (1992) 
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6200 6400 6600 6800 7000 7200 7400 

Wavelength (A) 



Fig. 2. — Continuum subtracted narrow-slit AAT spectra of SN 1987A and its ring structure showing the spectral region between 
6200 A and 7400 A on three epochs, Jan. 1991 (1416 days), Dec. 1994 (2864 days), and Dec. 1998 (4309 days). The scale of the plots has 
been selected to better show the faint lines measured from these narrow slit spectra whereas the brightest lines (measured from the wide 
slit spectra) extend outside the plotting range. The spectral lines included in this study have been identified. 



that there are indeed no forbidden lines present in the 
Star 3 spectrum. 

WFPC2 images obtained using the F502N and F656N 
filters were selected from the HST archive at six epochs 
contemporary with our AAT spectra (e.g., see Fig. 1). 
Pipeline reduced WFPC2 data were acquired from the 
HST archive and tasks from the IRAF STSDAS package 
were used to combine the pipeline reduced repeat im- 
ages rejecting for cosmic rays and hot pixels. Aperture 
photometry was then performed for Star 3 using the rec- 
ommended aperture diameter of 0.5" and then corrected 
to an infinite aperture with the standard aperture correc- 
tion. As we are interested in the Ha emission line flux on 
top of the stellar continuum, the narrow band fluxes of 
Star 3 needed to be corrected for the contribution from 
the continuum. The level of the continuum at the wave- 
length of Ha was measured from the STIS spectra and 
was subtracted from the measured F656N flux densities. 
However, as Star 3 is known to vary with an amplitude 
of ^0.5 magnitudes at optical wavelengths its continuum 
level can be expected to vary accordingly. Therefore, the 
continuum subtraction can be expected to introduce an 
error of about ±20% to the Ha emission line fluxes. The 
continuum subtracted Ha fluxes are listed in Table 2. As 
no forbidden emission lines have been seen in the Star 3 
spectrum the F502N band (centered on the [O III] line 



at 5007A) can be expected to contain pure continuum 
flux from the star. Therefore, the narrow band fluxes 
observed with the F502N filter can be used to estimate 
the B-band magnitudes of Star 3. The magnitudes of 
Star 3 transformed to the Johnson B-band are listed in 
Table 2. In the same way a magnitude of = 15.01 was 
obtained for Star 2 which is close to the earlier ground 
based measurement = 14.97 (Walborn et al. 1993). 

Average Ha and H/3 fluxes of 4.4 and 1.3 x 10~ 14 erg 
s _1 cm -2 , respectively, were found for Star 3, with a 
maximum variability of ^20% over a 6 year time span. 
These fluxes were used to correct for Star 3 contamina- 
tion of the hydrogen line fluxes in the wide slit AAT ob- 
servations. The additional uncertainty in the Ha and H/3 
fluxes due to the correction for the Star 3 contamination 
was estimated to be ±5% for the 1416 and 1680 days 
epochs and ±15% for the 2864 and 4309 days epochs. 
The absolute flux calibration for the vignetted narrow 
slit spectra was obtained by scaling them to match the 
unvignetted wide slit spectra using the brightest forbid- 
den lines present in both the spectra. Deriving line fluxes 
for the whole inner ring this way assumes homogeneous 
ring geometry Comparison to the line fluxes from the 
wide slit spectra indicate an uncertainty of up to ±40% 
for the line fluxes measured from the narrow slit spectra. 

The accuracy of the absolute fluxing of our AAT ob- 
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TABLE 3 

Line fluxes of the inner ring from the literature (307-2122 days). 



line 307 511 552 574 586 668 678 695 734/735 840 958 1050 1114 1280 1344 1348 1469 1734 1822 2122 



Ha 






... 283 


... 269 ... 


324 




251 




143 


140 












H/3 


17 


46 64 


... 72 


... 74 ... 


100 




69 




42 


39 












H7 


10 


... 27 


... 35 


... 24 ... 


21 




29 




14 


16 












H<5 


25 




... 26 


8 


14 








8.8 


9.1 












H(Pa 7 ) 












35 




50 






28 


10 


12.4 


9.5 


6.0 


H(Pa/3) 










50 


40 90 




70 






34 


10 


9.0 


7.0 


6.5 


H(Brj) 
























5.5 


1.5 


2.0 


1.5 


He II A4686 




... 24 


... 19 


... 14 ... 


24 




14 




5.0 


6.8 












He I A5876 






... 17 


... 19 ... 


15 




16 




7.1 


7.0 












He I A7065 


















1.8 














He I A10830 






1400 ... 


2400 ... 2200 


1500 


1000 1400 




650 






390 


360 


120 


90 


53 


He I A20581 












5 




5 








7.0 


3.5 


3.0 


2.0 


[N II] AA6548.6583 






... 792 


... 796 ... 


938 




842 




286 


349 












[N II] A5755 






20 


24 


23 




20 




9.5 


8.0 












[O I] AA6300.6364 










11 








3.6 














[O II] AA3726,3729 










77 




75 




71 


49 












[O II] AA7319J331 


















12.9 














[O III]AA4959,5007 


433 


... 973 


... 857 


... 668 ... 


825 




393 




150 


144 












[O III] A4363 


27 


24 35 


42 


24 






18 




5.5 


5.9 












[Ne III] A3869 










69 




60 




18 


19 












[S II] AA6716.6731 








... 36 ... 


41 




38 




17 


20 












[S II] AA4069,4076 






... 24 


... 15 ... 


15 








4.5 


4.1 












[S III] AA9069,9532 
























61 


40 






[Fe II] A5159 


















0.8 














[Fe II]A7155 


















0.8 














[Fe II]A12567 
























6.0 


4.5 


3.0 


2.0 


[Fe II] A16435 
























9.0 


4.0 


3.0 


1.5 


[Fe III]A4658 


















1.9 















The fluxes are in units of 10 15 erg s 1 cm 2 . The fluxes for 307 days were taken from Wampler & Richichi (1989), and were reddened 
back to observed fluxes using E(B — V) = 0.05 and 0.10 for the Galaxy and the LMC, respectively, and the extinction laws from Howarth 
(1983) assuming R v = 3.1. The fluxes for 511, 552, 586, 678, and 735 days are from Wampler et al. (1989). The fluxes at 574, 668, 
695, 734, 840, 958, and 1114 days are from Meikle et al. (1993). The fluxes at 1050 days are from Menzies (1991), reddened back to the 
observed fluxes using E(B — V) = 0.05 and 0.19 for the Galaxy and the LMC, respectively, and the extinction laws from Howarth (1983). 
The [O III] (4363, 4959, 5007) and [N II] (5755, 6548, 6583) fluxes given by Menzies (1991) are from day 1050 whereas the other line fluxes 
are averages of the observed fluxes between ~ 850 and 1200 days. The fluxes at 1280 days are from Wang (1991). The fluxes at 1344 days 
are from Cumming (1994), reddened back to the observed fluxes using E(B — V) = 0.05 and 0.16 for Galaxy and the LMC, respectively 
and the extinction laws from Howarth (1983). The fluxes at 1348, 1469, 1734, 1822, and 2122 days are from Fassia et al. (2002). 



servations were also checked using the HST inner ring 
data (see also Mattila 2002). The F502N filter is very 
useful for this as it contains the well isolated [O III] line 
at 5007 A. In addition, we used the F656N imaging data 
as this filter covers well the Ha and the [N II] line at 
6548 A, and is only slightly contaminated by the nearby 
[N II] line at 6583 A. Therefore, the fluxes of the inner 
ring were measured from the F502N and F656N images 
using an aperture with a 1.35" radius, an ellipticity of 
0.75, and a PA selected to match the orientation of the 
ring in a given HST image. The flux from the SN in the 
middle of the ring was measured with a circular aperture 
with a radius of 0.32" selected to include most of the flux 
from the SN avoiding the emission from the ring. The 
flux measured from the SN was typically ^5-10% of the 
ring flux in the F656N images. The fluxes of the inner 
ring, corrected for the SN emission, are listed in Table 2 
together with the corresponding inner ring fluxes as mea- 
sured from the AAT spectra and interpolated (after the 
Star 3 correction) to the epochs of the HST observations. 
Comparison between these fluxes shows the AAT obser- 
vations to be well consistent with the HST measurements 
to within a maximum deviation of 25%. 



2.3. VLT Spectroscopy 

High resolution (FWHM ~ 6 km s _1 ) spectroscopic ob- 
servations were obtained with the Ultraviolet and Visual 
Echelle Spectrograph (UVES) at the ESO VLT during 
2000 December 9-14 (-5040 days post-explosion), 2002 
October 4-7 (~5703 days post-explosion), 2005 March 
21 - April 12 (~6618 days post-explosion), 2005 Octo- 
ber 20 - November 12 (^6826 days post-explosion), 2006 
October 1 - November 15 (~7183 days post-explosion), 
and 2007 October 23 - November 28 (~7559 days post- 
explosion). The UVES data cover the spectral range 
~3100-10,000 A. The slit width was 0.8" and oriented 
at a position angle of 30° across the SN. The airmass of 
these observations ranged between ~1.4 and 1.8 and the 
seeing spanned 0.4-1.4". The total exposure time was 
about 10,000 sec in each epoch. For full details of the 
observations see Table 1 in Groningsson et al. (2008b). 
An atmospheric dispersion corrector was always used for 
these UVES observations to minimize any wavelength 
dependent slit losses. The data reduction, flux calibra- 
tion and the measurement procedure of the line fluxes are 
discussed in detail in Groningsson et al. (2008a). With 
the relatively narrow slit only a fraction the ring flux 
was covered by the observations. To estimate this frac- 
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tion and to check the accuracy of the flux calibration we 
compared the UVES fluxes to data taken with the HST 
at a similar epoch. From this comparison we estimated 
that the correction factor needed to compensate for the 
slit losses was 2.4 (Groningsson et al. 2008b, Appendix 
A) and the uncertainty of the absolute fluxes should be 
less than —±20% (see also Groningsson et al. 2008a). 

2.4. The Line Light Curves 

We assembled from the literature all the available 
optical/near- IR line fluxes of the inner ring between 
-300 and -2100 days. Early time optical CSM line 
fluxes were collected from Wampler & Richichi (1989), 
Wamplcr et al. (1989), Wang (1991), Menzies (1991), 
and Cumming (1994). Near-IR (NIR) fluxes covering 
both early and late time epochs were taken from Mcikle 
et al. (1993), and Fassia et al. (2002). The fluxes from 
Wampler & Richichi (1989), Menzies (1991), and Cum- 
ming (1994) were reddened back to the observed fluxes 
to enable consistent dereddening for the data from all the 
epochs. The fluxes given by Cumming (1994) had been 
dereddened using the extinction laws of Howarth (1983) 
(R. Cumming, private communication). Without any 
additional information, the same extinction laws were 
also used when reddening the fluxes given by Wamplcr 
& Richichi (1989) and Menzies (1991). In all cases it was 
assumed that the Galaxy contributed E(B — V) = 0.05 
with the rest of the reddening (see the caption of Table 
3) due to the extinction within the LMC. Meikle et al. 
(1993) estimate errors of ±15% for their He I A10830 line 
fluxes. Also, Fassia et al. (2002) estimate errors of ±35% 
for their IJ-band line fluxes. However, the other authors 
do not give any error estimates for their measurements. 
We therefore adopted errors of ±30% for these early time 
optical line fluxes. The fluxes from the literature are 
listed in Table 3. 



The high resolution VLT spectra were used to check 
if any of the lines observed in the much lower resolu- 
tion AAT spectra suffered from contamination by an- 
other line. For some of the lines (e.g., [O II] AA7319,7331 
and [Ca II] A7324) the spectral resolution of the AAT 
narrow-slit observations was sufficient to allow deblend- 
ing the lines using Gaussian fitting. However, the He I 
A3889 and [Ne III] A3968 lines were found to be strongly 
blended with nearby lines of hydrogen and were therefore 
excluded from the analysis. The measured AAT and VLT 
optical line fluxes at 1416, 1680, 2864, 4309, 5040 days 
are listed for H, He, N, O, Ne, S, Ar, Ca and Fe in Ta- 
ble 4. The Ha and H/3 fluxes at 1416, 1680, 2864, 4309 
days were corrected for contamination by Star 3. The 
fluxes measured directly from the wide slit spectra were 
found to have uncertainties of ±15 - 21%, the fluxes mea- 
sured from the narrow slit spectra uncertainties of ±43%, 
and the fluxes measured from the VLT spectra uncer- 
tainties of ±20% except for [S III] AA9069,9532 in the 
red end of the spectrum for which ±50% was adopted. 
The correction of the Ha and H/3 fluxes from epochs 
earlier than 1400 days was not possible without know- 
ing the details of the observations and the subsequent 
data reductions. However, Star 3 contamination at these 
early epochs is not significant, since the ring was much 



brighter than Star 3. To investigate the effects of the col- 
lision between the SN ejecta and the inner ring material 
we also measured the fluxes for [O III] AA4959,5007, [Ne 
III] A3869, [Ne IV] AA4714,4726, [Ar III] AA7136.7751, 
[Ar IV] AA471 1,4740 and [Fe VII] A5721 at 5703, 6618, 
6826, 7183, and 7559 days from the VLT spectra. These 
fluxes are listed in Table 4. 

3. MODELING THE LINE LIGHT CURVES 

The latest version of the photoionization code used and 
described in Lundqvist & Fransson (1991), Lundqvist 
(1992, 1999, 2007) and LF96 was used here to model 
the emission line fluxes from the inner ring around SN 
1987A. In the model the inner ring is initially ionized 
by the soft X-ray and UV photons emitted in the SN 
shock break-out. We have used the 500fulll model of 
Ensman & Burrows (1992) for the ionization. It was 
confirmed by Blinnikov et al. (2000) to give a good rep- 
resentation of the break-out burst (cf. Lundqvist 2007). 
Initially the ring material was assumed to be optically 
thick to the ionizing flash and therefore, the outer radius 
of the ionized gas determined by the number of ionizing 
photons and the density of the gas, i.e., the emission re- 
gion in the ring was assumed to be ionization bounded. 
The outer radii of the different model density compo- 
nents were checked and if needed truncated to match the 
observed dimensions of the inner ring (see Sect. 3.1). 
The SN flash sets up the initial ionization structure of 
the ring, and the gas then recombines and cools. The 
model includes cooling by ions of H, He, C, N and O 
with all the ionization stages, and Ne, Na, Mg, Al, Si, 
S, Ca, Ar and Fe with the first 9, 10, 11, 12, 11, 11, 11, 
11 and 16 ionization stages, respectively. The code fol- 
lows the ionization, recombination and cooling of the gas 
with time including collisional and Auger ionization (i.e., 
a high energy photon ejecting an inner shell electron fol- 
lowed by ejection of an outer shell electron), radiative 
and dielectric recombination, collisional excitation and 
de-excitation, and charge transfer. As in LF96, the ge- 
ometry of the ring is modeled as a section of a spherical 
shell with the inner and outer edges at the same opening 
angle (8 — 2.5°) seen from the SN located in the middle 
of the circular ring. The ring is divided into 70 concentric 
"shells", the outer ones being more massive, (i.e., bigger 
volume) than the inner ones. This is because the ther- 
mal and ionization structure changes more rapidly in the 
inner parts of the ring thus requiring better "sampling" . 

The input parameters of the code are the mass and ra- 
dius of the ring, the gas number density, the spectral and 
temporal characteristics of the ionizing SN flash, and the 
elemental abundances in the ring. The program extracts 
emissivities for the most common UV-IR emission lines 
of H, He, C, N, O, Ne, Na, Mg, Al, Si, S, Ar, Ca and Fe at 
several ionization stages. For many of the ions, multilevel 
model atoms are used to calculate the emission. While 
the levels of ionization are treated time-dependently, the 
levels within the model atoms are in steady state, ex- 
cept for H and He where the population of all atomic 
states are treated time-dependently, with a total of 56 
levels for H I-II and 72 levels for He I-III. Strong reso- 
nance lines of all ions, and selected lines between excited 
states, are treated with escape probabilities, taking the 
limited thickness of the ring into account. In the multi- 
level atoms of H and He, photoionization is assumed to 
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TABLE 4 

Line fluxes of the inner ring from AAT and VLT (1416-7559 days) 



line 


1416 


1680 


2864 


4309 


5040 


5703 


6618 


6826 


7183 


7559 


Ha 


226'' 


186 b 


86 b 


76 b 


55 












H/3 
H7 




476 
14 b 


32 b 
14 


15 b 
8.4 


14 












H5 




7.9 




3.3 














He II A4686 




9.3 


5.7 b 


2.4 














He I A5876 


12 b 


15 


6.2 b 


4.0 


2.8 












He I A7065 


3.9 b 


4.5 


1.9 b 


1.0 6 


0.8 












[N II] AA6548.6583 


633 b 


580 b 


515 b 


340 b 


226 












[N II] A5755 


15 b 


15 b 


10 b 


7.3 b 


2.9 












[O I] AA6300.6364 


8.9 6 


15.1 


18 b 


14 b 


9.9 












[O II] AA3726,3729 




99 


70 b 


45 


36 












[O II] AA7319,7331 


22 


16 


9.8 


5.5 


3.4 












[O III] AA4959,5007 




189 b 


61 b 


48 


28 


40 


53 


60 


57 


69 


[O III] A4363 




5.6 b 


3.2 


2.1 


1.4 












[Ne III] A3869 




18 


11 


12 




5.2 


8.5 


8.5 


7.8 


6.9 


[Ne IV] AA4714,4726 .. 


... 


... 


... 


... 


0.11 




0.10 


0.16 


0.20 




[S II] AA6716.6731 


40'' 


36 b 


49 b 


38 b 


28 












[S II] AA4069,4076 




5.6 


9.2 


6.0 


4.3 












[S III] AA9069.9532 


51 b 






14 


5.0 












[Ar III] AA7136J751 


7.0 


7.7 






1.2 












[Ar IV] AA4711,4740 












1.3 


1.7 


1.9 


1.8 


1.6 


[Ca II] AA7291J324 


5.7 


6.8 


8.0 


6.2 


3.3 












[Fc II] A5159 




0.9 


1.0 


1.1 


0.48 












[Fc II] A7155 


1.2 


1.5 


1.9 


1.0 


0.73 












[Fc II] A7453 


0.6 


0.7 


0.7 


0.2 


0.24 












[Fe III] A4658 




2.2 


1.1 


1.2 














[Fe VII] A5721 












0.14 


0.10 


0.092 


0.11 


0.096 



The fluxes are in units of 10~ 15 erg s _1 enr 2 and are all observed fluxes (not corrected for reddening) from this study. The fluxes at 
1416-4309 days are from the AAT and the later ones from the VLT. The H« and H/3 fluxes for these epochs have been corrected for the 
contamination by Star 3. The fluxes measured directly from the AAT wide slit spectra (marked with b) have a smaller uncertainty (±15 - 
21%) than the fluxes measured from the narrow slit spectra (±43%). The fluxes at 5040-7559 days have uncertainties of ±20% except [S 
II] AA9069,9532 for which the uncertainty is ±50%. 



TABLE 5 

Atomic data for Fe II, Fe III and Fe VII. 



Ion 


Levels Temperature 


((X1000K) A value ref. Coll. strength rcf. 


Fe II 


46 1-20 


Nussbaumer & Storey (1988) Bautista & Pradhan (1996) 
Garstang (1962) Bautista, private communication 


Fe III 


27 3-70 


Nahar & Pradhan (1996) Zhang (1996) 


Fc VII 


9 2-1000 


Berrington et al. (2000); Young ct al. (2005) Bcrrington ct al. (2000) 



occur only from the ground state. While this is a good 
approximation to calculate the degree of ionization, the 
population of some excited states with low de-excitation 
rates could be less well modeled. 

Although the ionization fractions for the 16 first ioniza- 
tion stages of iron are calculated in the photoionization 
code, the code does not calculate individual line emissiv- 
ities for many of the iron ions. Instead, the cooling was 
pre-calculated using multilevel atoms and tabulated for 
various densities and temperatures as described in LF96. 
The atomic data included in the code are not sufficiently 
accurate to calculate the populations of individual energy 
levels with an accuracy allowing the iron line emissivi- 
ties to be determined. Therefore, we stored temperature, 
electron density and levels of ionization for iron at sev- 
eral time steps and ran a separate code to determine the 
line fluxes. For this we collected atomic data for Fe II, Fe 
III and Fe VII from the literature (see Table 5), to allow 
accurate modeling of the [Fe II] , [Fe III] and [Fe VII] line 



light curves. Thanks to the Iron Project (Hummer et 
al. 1993), these atomic data have improved significantly 
and make such abundance determinations more reliable. 
For Fe II, a 46-level model atom was used. We took the 
Einstein A-values for IR transitions from Nussbaumer & 
Storey (1988) and the other transitions from Garstang 
(1962). The more recent A-values by Quinet et al. (1996) 
were not used since they have been observed to yield dis- 
crepancies when used in modeling (Bautista & Pradhan 
1998). The collision strengths were taken from Bautista 
& Pradhan (1996), and Bautista, private communication. 
As the temperature varies strongly with time and across 
the ring, 8 different sets of collision strengths are used 
here for temperatures ranging between 1000 and 20,000 
K. For Fe III, a 27-level model atom was used. The Ein- 
stein A values were taken from Nahar & Pradhan (1996), 
and the collision strengths from Zhang (1996). Here 16 
different sets of collision strengths were used for temper- 
atures ranging from 3000 to 70,000 K. For Fe VII we 
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used a 9-level model adopting data from Bcrrington et 
al. (2000) and the corrected A-values by Young et al. 
(2005). 

Finally, equipped with the level populations of the ions 
of interest the total line fluxes observed from the whole 
ring were obtained by integrating the line emissivities 
over the volume of the ring (see Lundqvist & Fransson 
1991 and LF96 for details). This takes into account the 
time delay caused by the physical size of the ring and the 
finite velocity of light. 

3.1. The Density Structure and Abundances 

The dimensions of the inner ring used for the model- 
ing were obtained from the observations by the HST. The 
thickness of the ring has been constrained by HST imag- 
ing (Jakobsen et al. 1991; Plait et al. 1995) to be around 
10-15% of the ring radius in the plane of the sky. This 
corresponds to ~20% of the ring radius for a circular ring. 
The ionized gas of the lxlO 3 cm" 3 density component 
extends out to 1.26 x the ring inner radius, which we 
have assumed to be 6.2 x 10 17 cm (cf. Lundqvist 1999), 
and therefore the outer radius of this density component 
was truncated at 7.4 x 10 17 cm to match the observed 
dimensions of the inner ring. The elemental abundances 
(by number of atoms) used for the modeling here were 
initially as follows: Hc/H = 0.25, N/C = 4.5, and N/O = 
1.0 with an overall mctallicity, Z, of 0.35 times solar, i.e. 
Z/(H+Hc+Z) = 0.35 times the solar ratio. The observed 
line fluxes were corrected for extinction using the same 
extinction laws for the LMC and the Galaxy as in LF96, 
assuming an Eb-v of 0.05 and 0.15 for Galaxy and the 
LMC, respectively. The extinction laws from Fitzpatrick 
(1985), and Savage & Mathis (1979) were used for the 
LMC and the Galaxy, respectively. The dereddened line 
fluxes were initially compared with the models by eye. 
The mass and density structure of the ionized gas in the 
inner ring was determined by searching for satisfactory 
model fits for (1) the absolute line fluxes of Ha and H/3, 
and (2) the observed time evolution of the line fluxes of 
all the other elements. Guided by the previous studies 
(LF96), 3xl0 4 atoms cm -3 was selected as the highest 
density component. Lower density components were also 
discussed in LF96, and the lowest density determined so 
far is that of Mattila (2002), Mattila et al. (2003) and 
Groningsson et al. (2008b) who find indications for den- 
sities down to 10 3 atoms cm -3 . We indeed find that 
at least two additional components, lxlO 3 and 3-5 xlO 3 
atoms cm" 3 , were required, in addition to that of 3x10 
atoms cm -3 , to produce the observed line light curves. 

To decide on the weights (i.e., the masses of ionized 
gas) for each density component we sought for the op- 
timal set of weights to give the lowest x 2 value for the 
model fit simultaneously to Ha, [N II] AA6548,6583, [O 
III] AA4959,5007 and [S II] AA6716,6731. The line fluxes 
at 1280 and 1344 days appear to be systematically lower 
than the fluxes from the other authors at similar epochs 
and therefore we decided not to include these points in 
the x 2 fits. Also, we were not able to find a satisfactory 
fit for the H/3 line fluxes at 307 days and decided to con- 
centrate on the epochs later than this for all the lines. 
We found that the optimal fit was obtained using 0.65, 
1.30 and 0.45 as the weights for the 1 x 10 3 , 3 x 10 3 
and 3 x 10 4 atoms cm" 3 density components, respec- 



tively. A slightly higher x 2 value was found for the 1 x 
10 3 , 5 x 10 3 and 3 x 10 4 atoms cm" 3 combination of 
densities. The masses of ionized gas in the three density 
components were about 1.2, 3.4, and 1.2 x 10 " 2 M Q , re- 
spectively. The lowest-density components dominate the 
line light curves at the later times. High resolution HST 
images indicate that this low-density gas is situated, on 
average, further away from the SN than the gas with a 
higher density (Lundqvist & Sonncborn 1997). However, 
in reality there is probably a continuous range of den- 
sities within the ring rather than a few discrete density 
components. 

The elemental abundances yielding the observed ab- 
solute line fluxes were then determined. This was done 
by adjusting the total number of emitting ions relative 
to hydrogen. The effects of the modified abundances 
on the temperature and ionization structure of the ring 
were also included by rerunning the photoionization code 
after initially finding the required abundances. The op- 
timal abundances were then determined by searching for 
the number of emitting ions yielding the lowest x 2 value 
for each line (see Table 6). The resulting model fits are 
shown in Figs. 3-6 with the contributions of the different 
density components indicated. 

3.2. UV Lines and Average Abundances 

The UV emission line fluxes from the inner ring were 
modeled by LF96 up to ^2000 days from the explosion. 
They found gas densities of 6 x 10 3 , 1.4 x 10 4 , and 3.3 
x 10 4 atoms cm" 3 , and derived abundances for He, C, 
N and O (see Table 7). In Fig. 7 we compare the model 
fluxes with the observations of a few UV lines to test 
how well our selected model parameters can explain the 
UV data. For these plots the same density components 
were used as for the optical lines whereas the abundances 
were adopted from LF96. We find that for the UV lines of 
helium and oxygen our model parameters and the abun- 
dance estimated in LF96 give reasonable fits, especially 
at the later epochs, giving more confidence in our abun- 
dance estimates for these elements. However, to satis- 
factorily model the UV lines of carbon and nitrogen a 
different combination of densities and abundances would 
be required. For both the elements the line light curves 
were successfully modeled in LF96 and in Lundqvist et 
al. (2007). However, for modelling the N V A1240 line 
light curve LF96 found that resonance scattering in the 
circumstellar gas external to the inner ring was required. 
This is not included in the models used for this study. 
We note that the different abundances derived in LF96 
and here should provide a useful estimate of the uncer- 
tainty in the abundances. We therefore decided to also 
include the abundances based on all the UV lines studied 
in LF96 except N V A1240 when estimating the average 
abundances with their associated uncertainties. There- 
fore these abundances are also included in Table 6 to- 
gether with the estimates based on the optical lines. 

Based on the data in Table 6, the average abundances 
and their standard deviations were estimated for each 
element and these are listed in Table 7. However, we de- 
cided to exclude a few lines when calculating the average 
abundances. These were He I A10830, [O I] AA6300,6364, 
and the near-IR [Fe II] line at 1.6435^m. The He I 
A10830 data yielded He/H = 0.54 which is -3 times 
higher than obtained from the other He I lines. We note 
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Fig. 3. — The observed fluxes of lines of H and He plotted together with the best fit (see Table 6) model fluxes. The observed line fluxes 
have been dereddened assuming Eg_y of 0.05 and 0.15 for the Galaxy and the LMC, respectively. The long dashed, short dashed and 
dotted lines correspond to the density components of 3xl0 4 , 3xl0 3 , and lxlO 3 atoms cm -3 , respectively. The sum of these models is 
shown as a solid line. 



that a similar very high helium abundance of ^0.45 was 
also obtained by Allen, Meikle & Spyromilio (1989) from 
their narrow line He I A10830 observations at 400-700 
days. Fig. 3 shows that He I A10830 flux is underpro- 
duced in our models at epochs earlier than 1500 days, but 
overproduced later being more consistent with the abun- 
dances derived from the optical He I lines. Including only 
the three last points in the fitting yields He/H = 0.23 
which is similar to the estimates obtained from the other 
lines. This could point at problems with fluxing in the 
earlier He I A10830 data. It is, however, well-known that 
He I A10830 is notoriously difficult to model accurately 
because of the highly metastable 2 3 S level which allows 
for optical depth effects to become important for transi- 
tions to this level (e.g., Osterbrock & Ferland 2008). We 
have included this in our models and are able to boost 



the He I A10830 / He I A5876 ratio by 10 - 20% com- 
pared to the numbers in Porter et al. (2005). The line in 
our sample which is most attenuated by optical depth ef- 
fects is He I A7065, and the fact that our models give the 
same helium abundance for the two He I AA5876, 7065 
lines suggests that optical depth effects are included ad- 
equately. We caution that photoionization from excited 
states (cf. above), is not included in our models, but it 
is not clear how that could boost He I A 10830 and not 
other transitions to the 2 3 S level. We, however, leave 
the more detailed investigation of the He I A10830 line 
light curve for a future study and do not include it here 
in our estimate for the helium abundance. The near-IR 
[Fc II] line at 1.6435^m is also excluded due to fluxing 
difficulties at these wavelengths (see Fassia et al. 2002). 
For the [O I] AA6300,6364 line, we found that our com- 
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Fig. 4.— The observed fluxes of [N II], [O I], [O II] and [O III] lines plotted together with the best fit (see Table 6) model fluxes. The 
observed line fluxes have been dereddened assuming Eg_y of 0.05 and 0.15 for the Galaxy and the LMC, respectively. The long dashed, 
short dashed and dotted lines correspond to the density components of 3xl0 4 , 3xl0 3 , and lxlO 3 atoms cm -3 , respectively. The sum of 
these models is shown as a solid line. 



bination of densities could not be used to obtain a satis- 
factory model fit. Furthermore, our fit to this line yielded 
an oxygen abundance ~4 times (or 15a) higher than the 
average abundance found using five lines of O II and 
O III. Assuming that the ring contains high-density re- 
gions that are surrounded by low-density gas, some of 
this low-density gas will be shielded from the most in- 
tense radiation from the SN flash and will therefore con- 
tain more neutral, or at least less ionized gas than if no 
shielding was taking place. In fact, since mid-1990 HST 
images have shown many " hot spots" appearing in differ- 
ent locations around the entire inner ring (e.g., Lawrence 
et al. 2000; Sugerman et al. 2002) as a result of the SN 
ejecta impacting "fingers" of dense ring material protrud- 
ing towards the SN (see Pun et al. 2002). In our model 
shielding by such higher density regions is not included 



and all the density components see the SN flash the same 
way. We expect the higher ionization lines of oxygen to 
be less affected by the effects of shielding and therefore 
yield better model fits and more reliable abundance es- 
timates in this study. The shielding effect should also 
affect other low- ionization lines (cf. below). 

3.3. Late Epochs 

As discussed in Lundqvist (1999, see also Chevalier & 
Dwarkadas 1995), gas in the H II- region interior to the 
inner ring, but probably also extending to larger radii 
than the ring at higher latitudes, can contribute to the 
observed line emission at late epochs. The density of 
this gas is most likely not higher than ~ 1.5 x 10 2 atoms 
cm~ 3 . Such low densities can retain an appreciable frac- 
tion of highly ionized gas like N V even at epochs as late 
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Fig. 5.— The observed fluxes of [O III], [S II], [S III], [Ne III] and Ar [III] lines plotted together with the best fit (see Table 6) model 
fluxes. The observed line fluxes have been dereddcned assuming Eg_y of 0.05 and 0.15 for the Galaxy and the LMC, respectively. The 
long dashed, short dashed and dotted lines correspond to the density components of 3xl0 4 , 3xl0 3 , and lxlO 3 atoms cm -3 , respectively. 
The sum of these models is shown as a solid line. 



as 6000 days. The main uncertainty is how much of the 
H II-region has been swept up by the SN blast wave at 
late epochs. From Fig. 3 of Lundqvist (1999) we note 
that after 4000-5000 days this effect can introduce an un- 
certainty larger than a factor of two for the modeled line 
emission from the H II-region. As opposed to the inner 
ring, the H II-region was already partially ionized prior 
to the SN explosion by the progenitor star. It therefore 
will not be heated as efficiently by the supernova flash as 
the inner ring. Nevertheless, we have chosen to include 
a 10 atoms cm~ 3 component in our models in the same 
way as the higher density components i.e. we have not 
included the effects of the earlier ionization, which only 
affects ions with the ionization potential of H I, or lower. 
The supernova can ionize gas with a 10 2 atoms cm~ 3 den- 
sity out to 2.23 x the ring inner radius. We, truncated 



the outer radius of this component in the same fashion as 
for the 10 3 atoms cm -3 component (see Sect. 3.1) but to 
correspond to 30% of the ring inner radius, to agree with 
the discussion in Lundqvist & Sonneborn (1997). We 
compared the post-5000 day light curves with the models 
adding the 10 2 atoms cm~ 3 component to the combina- 
tion of density components found in Sect. 3.1 (see Fig. 
8). For [O III], [Ne III], [Ne IV], [Ar III] and [Ar IV] 
the average abundances found when fitting the pre-5000 
day light curves were adopted (see Table 7). However, 
for the [Fe VII] A5721 line light curve an iron abundance 
higher than the average by one standard deviation was 
required to satisfactorily model the data. After the trun- 
cation the 10 2 atoms cm -3 component was multiplied by 
6 (corresponding to a mass of ionized gas of 1.8 x 10~ 2 
M to yield satisfactory fits to the observed line fluxes 
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Fig. 6. — The observed fluxes of [Ca II], [Fe II] and [Fe III] lines plotted together with the best fit (see Table 6) model fluxes. The 
observed line fluxes have been dereddened assuming Eg_y of 0.05 and 0.15 for the Galaxy and the LMC, respectively. The long dashed, 
short dashed and dotted lines correspond to the density components of 3xl0 4 , 3xl0 3 , and lxlO 3 atoms cm -3 , respectively. The sum of 
these models is shown as a solid line. 



at ^5000 days. While the contribution of the 10 2 atoms 
cm~ 3 component was negligible for the [O III], [Ne III] 
and [Ar III] lines, it was required to satisfactorily model 
the [Nc IV], [Ar IV] and [Fe VII] line light curves. As- 
suming a constant density between the inner ring radius 
and out to 1.3 x this radius, the low-density gas would fill 
a region with an opening angle of 7.2° above and below 
the ring plane in addition to the gas in the ring. 

At epochs later than ~ 5000 days a strong contribu- 
tor to the narrow line emission is the re-ionization of the 
ring by soft X-rays from the hydrodynamic interaction 
of the ring with the supernova ejecta. Groningsson et 
al. (2008b) demonstrate clearly how especially [O III] 
and [Ne III] lines increase in flux between ~ 5000 and 
6800 days after the explosion due to this effect. Later, at 
least some of these lines appear to level off, or decrease 



in flux. We have not attempted to model the hydro- 
dynamic interaction. However, we do, as suggested in 
Lundqvist (1999), from day 1150 include modest ioniza- 
tion to simulate the early part of the light curve of soft 
X-rays (Hasinger et al. 1996). We have extrapolated the 
linear increase of X-ray emission assumed in Lundqvist 
(1999) to late times. As shown in Park et al. (2007), 
the soft X-ray emission starts to increase much faster 
after ~ 4700 days than we have used in our modeling. 
We have refrained from trying to fine tune our models 
with detailed fits to the X-ray emission. This means 
that our models should underproduce the emission of at 
least some lines after ~ 5000 days. The amount of under- 
production is an indication of how much emission is the 
result of a shock precursor of the ejecta/blob collisions. 
In Fig. 8 this effect can be clearly seen for most of the 
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Epoch 

Fig. 7.— The observed fluxes of the He II A1640, C III] A1909, O III] A1665, N III] A1751, N IV] A1486, and N V A1240 UV lines plotted 
together with the model fluxes. The observed line fluxes (from Sonneborn et al. 1997) have been dereddened assuming Eg_^ of 0.05 and 
0.15 for the Galaxy and the LMC, respectively. The long dashed, short dashed and dotted lines correspond to the density components of 
3xl0 4 , 3xl0 3 , and lxlO 3 atoms cm , respectively, using the abundances found in LF96 (see Table 7). The sum of these models is shown 
as a solid line. 



lines considered here. 

4. DISCUSSION 

4.1. Abundances 

The main aim of this study was to determine the el- 
emental abundances of the gas in the inner ring of SN 
1987A. The results are summarised in Table 7, where 
we also list abundances found from fits to X-ray grat- 
ing spectra of the ring collision (Dewey et al. 2008, 
Zhekov et al. 2009), as well as for the LMC (Russell 
& Dopita 1992, Hughes et al. 1998, Hunter et al. 2007) 
and for solar elemental abundances (Anders & Grevesse 
1989, Grevesse et al. 2007). For the inner ring, LF96 
derived the abundances 11.40, 7.51, 8.26, 8.20 for He, 
C, N, O, respectively, in the same units as in Table 7, 
i.e., 12 + log[n(X)/n(H)] where n(X) is the number den- 



sity of element X. Those numbers were obtained from 
fits to mainly UV lines observed with IU E. In the cur- 
rent paper we find 11.20, 8.41, 8.27, for He, N and O, 
respectively. No carbon line exists in this optical data 
set. 

The average He/H-ratio found using the three optical 
helium lines considered in this study is 0.14 ± 0.04, and 
0.17 ± 0.06 when including also the He II A1640 UV line 
from LF96. The large error is due to the He II A4686 
line yielding ^2 x lower abundance than the rest of these 
lines. This helium abundance is lower than, but within 
the errors of, the 0.25 ±0.05 derived by LF96. We believe 
that our result is more robust than that of LF96 since 
their abundance was obtained from a single, excited He II 
line in the UV where the uncertainty in the reddening 
correction has a larger effect than in the optical. Hydro- 
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Fig. 8.— The observed fluxes of the [O III] AA4959.5007, [Ne III] AA3869.3968, [Nc IV] AA4714.4726, [Ar III] AA7136,7751, [Ar IV] 
AA4711,4740, and [Fe VII] A5721 lines plotted together with the model fluxes. For oxygen, neon and argon the abundances found when 
fitting the prc-5000 day light curves were adopted (see Table 7), and for iron an abundance higher than the average abundance by one 
standard deviation was used. The observed line fluxes have been dereddened assuming E^-v of 0.05 and 0.15 for the Galaxy and the 
LMC, respectively. The long dashed, short dashed, dotted and long-short dashed lines correspond to the density components of 3xl0 4 , 
3xl0 3 , 1x10 s and lxlO 2 atoms cm 3 , respectively. The sum of these models is shown as a solid line. 



gen and helium have also been updated in our models (cf. 
above) since LF96. He/H = 0.17 is <~ 2.0 times larger 
than the solar and the LMC values. Abundances derived 
from X-ray models of the ring collision (e.g., Zhekov et al. 
2009) have used the He and C abundances of LF96 and 
then varied other abundances to make best fits. With 
He/H = 0.17, instead of the 0.25 in LF96, the X-ray 
based abundances would need to be somewhat revised. 

The nitrogen abundance we find using only the opti- 
cal [N II] lines from this study is twice as high as the 
one found in LF96 whereas the oxygen abundances from 
these two studies are almost identical. The N/O-ratio 
based on only the optical lines is 1.8 ± 0.4, which is 
therefore significantly higher than the one found by LF96 
of N/O = 1.1 ± 0.4. Including also the UV lines from 



LF96 in our abundance estimates the N/O-ratio becomes 
1.5 ± 0.7, which is higher than, but still within the er- 
rors of, the one found by LF96 and the ratio N/O ~ 0.9 
found by the X-ray studies (Dewey et al. 2008, Zhekov 
et al. 2009). In Lundqvist (2007) we redid the calcu- 
lations from LF96 with the most recent version of the 
flash-ionization code, and obtained good fits to both the 
UV line light curves and the optical [N II] AA6548,6583 
and [O III] AA4959,5007 line light curves for epochs ear- 
lier than 1400 days for the same abundances as in LF96. 
The chief differences between the LF96 models and those 
in the current paper are mainly due to the difference in 
density distribution and that we now concentrate on data 
at later epochs. In LF96 we could not constrain gas with 
densities lower than 6 x 10 3 cm~ 3 , whereas our study here 
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TABLE 6 

Abundances for different lines 
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Best fit reduced x 2 values and the resulting abundances (relative 
to hydrogen) for different lines (see Figs. 3-6). Whether the line 
is included in the average abundance (see Table 8) is indicated in 
col. 4. 



is mainly sensitive to lower densities. We note that both 
the study of LF96 and X-ray studies are likely to probe 
the gas on the inside of the ring where the ring densities 
probably are higher than further out (cf. the discussion 
in Lundqvist & Sonneborn 1997). Here, we have con- 
centrated on the low-density gas and have assumed that 
this is not shielded from the UV flash by high-density 
gas closer to the supernova. This may not be true for 
all low-density gas and thus there could be a systematic 
difference between LF96 and the present analysis. Our 
estimated N/O-ratio is much higher than the range 0.03 
< N/O < 1 found for a sample LMC B-type supergiants 
by Hunter et al. (2007). However, a very high N/O-ratio 
of ^2.1 has also been found for the circumstellar nebula 
around the Galactic blue supergiant Shcr 25 by Hendry 
et al. (2008). Sher 25 has a spectral type similar to 
the progenitor of SN 1987A at the time of its explosion. 
Furthermore, Sher 25 is also surrounded by an hourglass 
shaped circumstellar nebula that is believed to have been 
ejected by the star about 6600 years ago and has been 
suggested to have many similarities to the ring structure 
of SN 1987A (Brandner et al. 1997). In the case of Sher 
25 Hendry et al. found also the photospheric N/O-ratio 
to be very high and consistent with the one they derived 
for the circumstellar nebula. 

Adopting the carbon abundance from LF96, together 
with Hc/H = 0.17 ± 0.06 and the N and O abun- 



dances found using both optical and UV lines in 
this paper, (N+0)/H = (4.7 ± 1.2) x 10~ 4 , the total 
(C+N+0)/(H+He) abundance is (4.3±1.0) x 10~ 4 . This 
can be compared with the old solar value of Anders & 
Grevesse (1989), (1.3 ± 0.1) x 10~ 3 , the new solar value 
of Grevesse et al., (7.1 ± 0.8) x 10~ 4 , and the LMC value 
according to Russell & Dopita (1992), (3.3 ±0.8) x 10~ 4 . 
In the recent paper of Hunter et al. (2007), the present- 
day LMC abundances of (C+N+0)/(H+Hc) for B-type 
stars is argued to be ~ 2.7 x 10~ 4 , i.e., consistent with 
the findings of Russell & Dopita (1992). Our results are 
roughly 1.6 times higher than found for the LMC, and 
about 60% of the modern solar values. The explanation 
could be due to the high helium abundance in SN 1987A. 
A simple test is to use LMC abundances of metals and 
helium, and to elevate these compared to hydrogen with 
the same factor as the Hc/H ratio in the inner ring of SN 
1987A, i.e., with a factor of ~ 2 to simulate the loss of hy- 
drogen from the surface of the progenitor. Doing so, the 
(C+N+0)/(H+Hc) abundance would be - 4.2 x 10~ 4 , 
i.e., nearly the same as we found for the inner ring of SN 
1987A. 

The iron abundance was estimated based on three 
[Fe II] lines and one [Fe III] line yielding 0.20 ± 0.11 times 
the solar value of Anders & Grevese (1989), which is con- 
sistent with the findings from the X-ray studies (Dewey 
et al. 2008, Zhekov et al. 2009) of - 0.2 x solar. It is also 
well within the range generally observed in the LMC, 
which is - (0.2 - 0.3)xsolar (Russell & Dopita 1992; 
Hughes et al. 1998). An iron abundance significantly 
lower than this could have indicated depletion onto dust 
grains in the ring material. Also, a silicon abundance 
of the X-ray emitting gas was found by Zhehov et al. 
(2009) to have a value typical for the LMC. However, 
recent Spitzer observations (Bouchet et al. 2006; Dwek 
et al. 2008) have found evidence for a significant amount 
of ~T0~ 6 M Q of silicate grains within the circumstellar 
ring. 

Our estimates for the neon, argon and calcium abun- 
dances are all based only on a single (doublet in the case 
of Ar and Ne) emission line and therefore, these abun- 
dances could be considered less reliable. However, for 
neon and argon we also compared late time (post 5000 
day) light curves of [Ne IV] A4714+4726 and [Ar IV] 
A4711+4740 lines with the models (see Fig. 8). For this 
their abundances found earlier (Table 7) were used and 
an additional 10 2 atoms cm~ 3 density component was 
included (see Sect. 3.2) giving a satisfactory match at 
^5000 days for all the six lines modeled. Also, the evo- 
lution of the [Ne III] doublet is expected to behave simi- 
lar to the strong and well-sampled [O III] doublet. This 
gives us more confidence on the estimated neon and ar- 
gon abundances and we conclude that they are probably 
accurate to a similar level as found for iron, i.e., ±~50%. 
Within these uncertainties we find that the abundances 
of both argon and neon are consistent with their LMC 
values and also with the abundances found in the X-ray 
studies of the inner ring. However, our abundances for 
calcium and sulphur are both higher than the ones found 
for the LMC and for the ring in the X-ray studies. We 
note that both the [S II] and [Ca II] lines used in our 
study are prone to the shielding effect discussed in Sect. 
3.2 for [O I], and we should thus only consider our abun- 
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TABLE 7 

Elemental abundances 



Element solar(AG89) solar(GAS07) LMC(+) LMC(X-ray,H) LMC(H07) 87A(X-ray,D) 87A(X-ray,Z) 87A(LF96) This work Lines 



He 


10.93(0.01,*) 


10.93(0.01,*) 
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Ar 


6.56(0.10) 


6.18(0.08) 


6.29(0.25) 








6.29 




6.23 


1 


Ca 


6.36(0.02) 


6.31(0.04) 


5.89(0.16) 








5.89 




6.51 


1 


Fe 


7.67(0.03) 


7.45(0.05) 


7.23(0.14) 


7.01(0.11) 




6.97(0.04) 


6.97(0.02) 




6 98+ - 19 


5 



The elemental abundances are in units of 12 + log[n(X)/n(H)] and the errors are given in parentheses. The abundance estimates from 
this work are the average of the estimates from different emission lines (see Table 6) and the quoted errors are the standard deviations. 
The number of emission lines used for each average abundance is given in col. 11. (*) = GAS07 = Grevesse, Asplund & Sauval (2007) 
(values are for photosphcric abundances). AG89 = Anders & Grevesse (1989). (+) = Russell & Dopita (1992). H = Hughes, Hayashi, 
Koyama (1998). D = Dewey et al. (2008). H07 = Hunter et al. (2007). Z = Zhekov ct al. (2009). LF96 = Lundqvist & Fransson (1996). 



dances for S and Ca as upper limits. 

4.2. Density Components and Late Epochs 

We found that satisfactory model fits could be pro- 
duced for the observed line light curves using a combina- 
tion of at least three different components with atomic 
densities of 1 x 10 3 , ^3 x 10 3 , 3 x 10 4 atoms cm~ 3 . 
The outer radius of the 1 x 10 3 component needed to 
be truncated to correspond to the observed dimensions 
of the ring whereas the higher density components were 
ionization bounded. The masses of ionized gas in the 
three density components were about 1.2, 3.4, and 1.2 
x 1CP 2 Mq, respectively, giving a total mass of ionized 
gas of ^5.8 x 1CP 2 M Q . For comparison, LF96 modeled 
the ring using three components 6 x 10 3 , 1.4 x 10 4 , and 
3.3 x 10 4 atoms cm -3 with masses of 2.4, 1.0, and 1.2 x 
10~ 2 M , respectively. Therefore, our total mass of ion- 
ized gas is ~ 25% higher than found by LF96, which is 
not surprising since we are here able to probe a broader 
range of density components than in LF96. 

We also compared late time (~5000 - 7500 days) light 
curves of [O III], [Ne III], [Ne IV], [Ar III], [Ar IV], and 
[Fe II] to the photoionization models and found that an 
additional 10 2 atoms cm -3 component of ~1.8 x 10~ 2 
M Q was required to explain the data at ^5000 days. Such 
low density gas is expected in the H Il-region interior to 
the inner ring which likely extends also to larger radii at 
higher latitudes. Assuming a constant density between 
the inner ring radius and out to 1.3 x this radius, the 
low density gas would fill a region with an opening angle 
of 7.2° above and below the ring plane in addition to 
the gas in the ring. The exact location and extent of 
this gas is not well constrained in our study, but we note 
that a thick H Il-region surrounding the equatorial ring 
as used in the models of Borkowski et al. (1997) can at 
least qualitatively explain the shape of the reverse shock 
of the ejecta/ring interaction (Michael et al. 2003). The 
H Il-region used in their model has an extent of ±30°, 
which is substantially larger than our ±7.2°. Assuming 
an opening angle of ~10° above and below the equatorial 
plane Zhekov et al. (2010) also found a density of ^10 2 
atoms cm~ 3 for the H Il-region based on their X-ray 
spectra. There could of course be gas with lower density 
than 10 2 atoms cm -3 at higher latitudes than in our 



model. Such gas would probably not contribute much to 
the strength of the narrow lines, but may still give rise 
to a similar structure of the reverse shock as inferred by 
Michael et al. (2003, see also Heng et al. 2006). At 
epochs later than ^5000 days our models underproduce 
the emission of most of these lines as expected due to the 
contribution from the interaction of the SN ejecta with 
the ring. 

5. CONCLUSIONS 

We have presented optical and near-IR line light curves 
of the inner circumstellar ring of SN 1987A covering the 
entire era from ^300 days when the ring first became 
visible at these wavelengths to ^5000 days when the col- 
lision between the SN ejecta and the CSM ring started 
to significantly affect these fluxes. We modeled the line 
fluxes with a photoionization code which follows the re- 
combination and cooling of the gas in the ring after the 
initial ionization by the SN flash. From this we draw the 
following conclusions. 

The absolute line fluxes of Ha and H/3, and the light 
curves for most of the lines studied here are modeled well 
with a combination of three density components lxlO 3 , 
^3xl0 3 and 3xl0 4 atoms cm~ 3 . The total mass of the 
ionized gas found was ^5.8 x 10~ 2 M Q with about 20%, 
60% and 20% of the mass in the three different com- 
ponents, respectively. Complementing our optical/near- 
IR abundances with estimates of Lundqvist & Fransson 
(1996) based on UV lines we found an He/H-ratio (by 
number of atoms) of 0.17 ± 0.06 and an N/O-ratio of 
1.5 ± 0.7. This helium abundance is roughly 30% lower 
than previously estimated but still ~2.0 times the so- 
lar and the LMC values. Our estimated N/O-ratio is 
higher than typically found for the B-type supergiants 
in the LMC. However, a very high N/O-ratio has also 
been found for Sher 25 (and its circumstellar nebula) 
which has sometimes been suggested to be a twin of the 
progenitor of SN 1987A. The total (C+N+0)/(H+He) 
abundance roughly 1.6 times its LMC value was obtained 
which is as expected assuming the abundances of metals 
to be enhanced by the same factor (^2.0) as He com- 
pared to their LMC values. Our iron abundance based 
on optical and near-IR lines is 0.20 ±0.11 times solar 
which is within the range of the estimates for the LMC 
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and therefore does not indicate a strong depletion onto 
dust grains in the ring material. 

An additional 10 2 atoms cm -3 component of ~1.8 x 
10 -2 M Q was required to explain the fluxes of a few high 
ionization lines studied at the late times. Such low den- 
sity gas is expected in the H Il-region interior to the inner 
ring which likely extends also to larger radii at higher lat- 
itudes. For a few lines we also presented light curves at 
epochs later than ^5000 days and found our models to 
underproduce the emission as expected due to the con- 
tribution from the interaction of the SN ejecta with the 
ring. The re-ionization of this gas by the X-rays from 
the interaction can therefore be expected to make the 
ring glow again especially in many high-ionization nar- 



row UV-lines making it an interesting target for future 
ground and space based observations. 
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